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Abstract 

In the context of a two-parameter (a, (3) deformation of the canonical commutation 
relation leading to nonzero minimal uncertainties in both position and momentum, 
the harmonic oscillator spectrum and eigenvectors are determined by using an exten- 
sion of the techniques of conventional super symmetric quantum mechanics combined 
with shape invariance under parameter scaling. The resulting supersymmetric part- 
ner Hamiltonians correspond to different masses and frequencies. The exponential 
spectrum is proved to reduce to a previously found quadratic spectrum whenever one 
of the parameters a, (3 vanishes, in which case shape invariance under parameter 
translation occurs. In the special case where a = (3 ^ 0, the oscillator Hamiltonian 
is shown to coincide with that of the g-deformed oscillator with q > 1 and its eigen- 
vectors are therefore n-q-boson states. In the general case where / a / /? ^ 0, the 
eigenvectors are constructed as linear combinations of n-g-boson states by resorting 
to a Bargmann representation of the latter and to (/-differential calculus. They are 
finally expressed in terms of a g-exponential and little g-Jacobi polynomials. 
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1 Introduction 



Studies on small distances in string theory and quantum gravity suggest the existence of a 
finite lower bound to the possible resolution of length Ax (see, e.g., [1, 2]) On the other 
hand, on large scales there is no notion of plane waves or momentum eigenvectors on generic 
curved spaces. It has therefore been suggested that there could also exist a finite lower 
bound to the possible resolution of momentum Ap (see, e.g., [3]). It is a natural, though 
nontrivial, assumption that minimal length and momentum should quantum theoretically 
be described as nonzero minimal uncertainties in position and momentum measurements. 

Such nonzero minimal uncertainties can be described in the framework of small correc- 
tions to the canonical commutation relation [4, 5] 



with a > 0, (3 > 0, and af3 < h 2 . In such a context, they are given by Ax = 



Since the canonical commutation relations lie in the very heart of quantum mechanics, 
studying the influence of small corrections to them in a quantum mechanical framework 
is interesting in its own right. It has indeed been argued [6] that such corrections may 
provide an effective description not only of strings but also of non-pointlike particles such 
as quasiparticles and various collective excitations in solids, or composite particles such as 
nucleons and nuclei. 

Solving quantum mechanical problems with the deformed canonical commutation rela- 
tion (1.1) may, however, be a difficult task. In the special case where a = and (3 > 0, the 
minimal uncertainty in the position turns out to be Ax = h\J~~j3~, whereas there is no nonzero 
minimal momentum uncertainty. As a consequence, equation (1.1) can be represented on 
momentum space wave functions (although not on position ones). Similarly, in the case 
where a > and (3 = 0, there is only a nonzero minimal uncertainty in the momentum 
Ap = and equation (1.1) can be represented on position space wave functions. On 
the contrary, in the general case where a > and (3 > 0, there is neither position nor mo- 
mentum representation, so that one has to resort to a generalized Fock space representation 
or, equivalently, to the corresponding Bargmann representation [4, 5, 7]. 

*In [4, 5], the two deforming parameters are denoted by a and (3. Here we reserve this notation for the 
dimensionless parameters to be introduced in section 2. 





K\J(3/(1 — h 2 af3) and Ap = hya/(l — h 2 aj3), respectively*. 
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The investigation of the harmonic oscillator with Hamiltonian 

H=^ + -mw 2 x 2 (1.2) 

where x and p satisfy the deformed canonical commutation relation (1.1), is an interesting 
and nontrivial topic. The eigenvalue problem for such an oscillator 

H\i, n ) = E n \i, n ) n = 0,1,2,... (1.3) 

has been solved exactly only in the case where a = and (3 > by using the momentum 
representation and the technique of differential equations [8]. This approach has been 
recently extended to D dimensions [9] and some ladder operators have been constructed [10]. 
While it is obvious that in the case where a > and (3 = 0, the eigenvalue problem can be 
treated in a similar way, that corresponding to both a > and (3 > is more complicated 
and, as far as we know, has not been solved so far. 

Over the years, it has been shown that supersymmetric quantum mechanics (SUSYQM) 
plays an important role in obtaining exact solutions of quantum mechanical problems (see, 
e.g., [11, 12]). In fact, all solvable quantum mechanical problems are either supersymmetric 
or can be made so. 

Among the various exactly solvable potentials, there is a certain class of potentials 
characterized by a property known as shape invariance [13]. Shape invariant potentials 
are potentials such that their SUSY partner has the same spatial dependence with possibly 
altered parameters. For such potentials, both the energy eigenvalues and the wave functions 
can be obtained algebraically without solving the differential equation [13, 14]. It turns 
out that the formalism of SUSYQM plus shape invariance (connected with translations of 
parameters) is intimately related to the factorization method developed by Schrodinger [15] 
and by Infeld and Hull [16] (for a comparison between these two methods and corresponding 
references see [12]). Other types of shape invariance, such as that connected with parameter 
scaling [17, 18], have provided a lot of new exactly solvable problems. 

The purpose of the present paper is to extend SUSYQM and the notion of shape invari- 
ance to eigenvalue problems in the context of the deformed canonical commutation relation 
(1.1) and to apply such an extension to the case of the harmonic oscillator with nonzero 
minimal uncertainties in both position and momentum (i.e., a > and (3 > 0). As a result, 
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we will derive exact results for the energy spectrum and the eigenstates of such a system 
in a purely algebraic way. 

This paper is organized as follows. The harmonic oscillator spectrum is obtained in 
section 2. Some special cases are reviewed in section 3. In section 4, the Hamiltonian 
eigenvectors are determined in explicit form. Section 5 contains the conclusion. Finally, 
some mathematical details are to be found in the appendix. 

2 Harmonic oscillator spectrum in the general case 

It is convenient to introduce dimensionless position and momentum operators, X = x/a 
and P = pa/h, where a = \Jhj (mwj is the oscillator characteristic length. They satisfy the 
commutation relation 

[X, P] = i(l + aX 2 + (3P 2 ) (2.1) 

where a and (3 denote the dimensionless parameters a = ah/{mui) and (3 = Pmhoj, respec- 
tively. Here we have a > 0, (3 > 0, and a(3 < 1. 

The dimensionless harmonic oscillator Hamiltonian is then given by 

h=^ = \{P 2 + X 2 ) (2.2) 

and the corresponding eigenvalue problem reads 

%n) = e„|Vn) n = 0,1,2,... (2.3) 

where e n = E n /(hw). To find the spectrum e n , n — 0, 1, 2, . . . , we shall proceed in two 
steps: we will show that the Hamiltonian h is factorizable, then we will prove that the 
factorized Hamiltonian satisfies a condition similar to the shape invariance condition of 
conventional SUSYQM. 

Let us first try to write h in the factorized form 

h = B + (g,s)B'(g,s) + e (2.4) 

where 

B ± (g,s) = ±={sX T igP) (2.5) 

and eo is the factorization energy. In (2.5), g and s are assumed to be two positive con- 
stants that are some functions of a, (3 and which go to 1 in the limit a, (3 — > 0. The 
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operators B + (g,s) and B~(g,s) are therefore counterparts of the standard harmonic os- 
cillator creation and annihilation operators, a + = (X — iP) / and a = (X + iP) / \/2, 
respectively. 

Inserting (2.5) in (2.4) leads to the equation 



h = \ 



(g 2 - (3gs)P 2 + {s 2 - ags)X 2 - gs + e . (2.6) 



For this expression of h to be equivalent to that given in (2.2), three conditions have to be 
fulfilled, namely 

g 2 -(5gs = 1 (2.7) 
s 2 -ags = 1 (2.8) 
e = \gs. (2.9) 

It can be easily shown that equations (2.7) and (2.8) admit positive solutions for g and 
s, given by 

g = sk s = - 1 = (2.10) 
yl — cxk 

where 

k = \{(3 - a) + y/i + l(j3- a )2. (2.11) 

In deriving the expression of s in (2.10), use has been made of the equivalence between 
the conditions 1 — ak > and aj3 < 1. Furthermore, one checks that k, g, s — > 1 for a, 
[3 — > 0, as it should be. We conclude that h can be written in the form (2.4), where the 
factorization energy eo is given in (2.9) and the two parameters a, f3 of the problem have 
been replaced by their combinations g, s, defined in (2.10) and (2.11). 
Let us now consider a hierarchy of Hamiltonians 

i 

h i = B + (g l ,s l )B-(g i ,s i ) + ^2^ i = 0,1,2,... (2.12) 

3=0 

whose first member coincides with (2.4), i.e., h = h. Here g iy Sj, 6j, i — 1, 2, . . . are assumed 
to be some positive parameters and g = g, s = s. We shall now proceed to prove that 
one can find values of g^ Sj, 6j, i — 1, 2, . . . , such that the condition 

B~(g i: Si)B + (g u = B + (g i+1 , s i+1 )B~(g i+u s i+1 ) + e i+1 (2.13) 
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is satisfied. Equation (2.13), written in operator form, is exactly equivalent to the equation 
used in the factorization method [15, 16] (see also [19]) and is similar to the shape invari- 
ance condition of conventional SUSYQM [11, 12, 13]. Note that the factorization method 
(and thus the SUSY one) is quite general (see, e.g., [19]) and can be used for finding the 
eigenvalues of arbitrary Hermitian operators with a bounded-from-below spectrum. As it 
is obvious that the eigenvalues of the Hamiltonian (2.2) are positive, we can apply the 
factorization method or the SUSY one with shape invariance to our problem. 
In explicit form, equation (2.13) reads 



(gf + (5g iSi )P 2 + (s 2 + ag iSi )X 2 + g iSi 



(9i+i- P9i+iSi+i)P 2 + (s 2 i+1 - ag i+1 s i+1 )X 2 - g i+1 s i+1 + e i+1 (2.14) 



where i — 0, 1, 2, This leads to the set of three relations 



9i+l ~ P9i+l s i+l 

$+1 ~ atg i+1 s i+1 



9i + P9iSi 
s 2 + agiSi 

\{9iSi + 9i+iSi+i). 



(2.15) 
(2.16) 
(2.17) 



At this stage, it is worth noting that by multiplying (2.15) by a and (2.16) by f3, then 
subtracting, we get the equation 



,,2„2 



9i+i ~ 7 s i+ i = 9i - 7 s 



2„2 



(2.18) 



where 
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(2.19) 



On iterating (2.18) and using (2.7) and (2.8), it is then obvious that the recursion for g iy 
Si runs along the arc of the hyperbola g 2 — 7 2 s 2 = 1 — 7 2 that lies in the first quadrant of 
the (g,s) plane. Since the equation of this hyperbola can also be written as uv — 1 — 7 2 , 
where u = g + 7s = and v = g — 7s = are the equations of the asymptotes, it may 
prove useful to replace (g, s) by (u,v). 

Let us therefore introduce the new combinations of parameters 



Ui 



9i + ISi 



9i ~ ISi 



(2.20) 
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The inverse transformation reads 

9i = \i.u i + v i ) s i = ^-(u i -v i ) (2.21) 
Z z7 

where the assumptions Qi, Si > impose that Ui > \v{\. On substituting (2.21) into (2.15), 
(2.16), then combining the two resulting equations, we obtain the relations 

u 2 +1 + qv 2 +1 = v 2 + qu 2 (2.22) 

u i+1 v i+ i = UiVi (2.23) 

where 

1 — \/ap 

and equation (2.23) coincides with (2.18). 

Equations (2.22) and (2.23) suggest a further transformation from Ui, Vi to 

di = u^i U = — (2.25) 

Ui 

where di and ti have the same sign, which is that of v^, and \ti\ < 1. According to (2.23), 
di is actually independent of i, 

d i = d = uv (2.26) 

which amounts to the hyperbola equation previously obtained. On the other hand, equation 
(2.22) can be rewritten as 

qk +l -U= q '^f ± ( 2 - 27 ) 
thus showing that qt i+i — ti = or 

U = q'H t = - = (2.28) 
u k + 7 

From (2.25), (2.26), and (2.28), we also obtain 

Ui = q i ' 2 u Vi = q~ i/2 v. (2.29) 

We conclude that the extended shape invariance condition (2.13) can indeed be satisfied 
by keeping the combination of parameters d constant while scaling the other combination 
of parameters t according to equation (2.28). Note that for % — > oo, Vi — > 0, so that the 
recursion actually reaches the hyperbola asymptote lying in the first quadrant of the (g, s) 
plane. 
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The eigenvalues e n of h are therefore given by 



n n— 1 | -y /n—1 | N 

Cn(q, t) = Y, e i = J2 9iSi + -g n S n = ~~ ^ ( U i _ U i ) + 9 ( M n ~ U n) 
i=0 i=0 Z 4 ^ \i=0 - / 

f « 2 ( [n] q + lq n ) ~ v 2 ( [n] s -i + V 



4 7 1 v l jy 2 
1 

where we have successively used equations (2.9), (2.17), (2.21), (2.29), (2.28), and the 
definitions 

M* = 7irr Wr^^3I = 9- n+1 W 9 - (2-31) 
In (2.30), we employ the notation e n (q, t) to stress that, apart from a multiplicative constant 
■u 2 /(47), the energy eigenvalues depend on the two parameters q and t, defined in (2.24) 
and (2.28), respectively. 

It can be easily shown that q(a, (3) = q(j3,a), u 2 (a, ff)/\4ry(a, ($)\ = u 2 (f3, a)/[47(/3, a)], 
and t(a, j3) = —t(/3, a), where we explicitly write down the dependence on the deformation 
parameters a, f3. As a result, the eigenvalues e n are symmetric under exchange of a and j3. 

From (2.30), we obtain for the ground state and excitation energies 

u 2 

e (q,t) = —(l-t 2 ) (2.32) 



and 

e.Aa. t) — e.nia. t) — 

2 V 1 n 

where 



e„fe,i)-eote,i) = iWl-4lM« n = 1, 2, . . . (2.33) 



* s "f£- (2M) 
Since g > 1 and t 2 < 1, the excitation energies grow exponentially to infinity when n — > oo. 

Such a feature has already been encountered before in SUSYQM and shape invariance 

associated with parameter scaling [17]. 

On using (2.5), (2.21), (2.28), and (2.29), the Hamiltonians (2.12) of the SUSYQM 

hierarchy can be written as 

hi = \{aiP 2 + hX 2 ) + a i = 0,1,2,... (2.35) 
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where 



u 2 , , / t 



(q l + t)[l + 



2(q+l)^ J \ g 1 " 1 



u 2 



h > = wwTT^-^y 1 -^) (2 - 36) 

For the supersymmetric partner of h — h , for instance, we obtain 

h = + (2/3 - a)k]P 2 + 1(1 + ak)X 2 + k) (2.37) 

1 — ak I z 1 2 J 

where fc is defined in (2.11). Going back to variables with dimensions, we get 

p 2 1 

Hi = huihi = h -miUJ 2 x 2 + Cihoj (2.38) 

2m,- 2 



with 

Tfi i 

rrii = — Ui = \JaibiU. (2.39) 

We conclude that the harmonic oscillator with nonzero minimal uncertainties in both 
position and momentum belongs to a hierarchy of Hamiltonians of the same type but 
with different masses and frequencies. The change of the kinetic energy term, appearing 
naturally in addition to the usual modification of the potential energy one, is a new feature 
arising from the deformation of the canonical commutation relation. It is distinct from the 
supersymmetric generation of combined potential and effective-mass variations that may 
be effected for Hamiltonians with both a position-dependent effective mass and a position- 
dependent potential in the context of conventional SUSYQM [20]. 

3 Some special cases 

In the present section, we will examine the two special cases where one of the parameters, 
e.g., a, vanishes or both parameters a, j3 are equal. 

3.1 Limit (\ — > 

We plan to show that in the limit a — > 0, the energy spectrum (2.30) reproduces the results 
previously obtained in the momentum representation [8, 9]. 
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For small a values, the parameter q behaves as 

q ~ 1 + 2^fa~f3 + 0(a) (3.1) 

so that 

q n ~ 1 + 2nya^ + 0(a) [n], ~ n + O(v^). (3.2) 

Furthermore, 

^(V-^l) ^ ffa +±/fc 2 (n-l) + 0(V5) 
i-^V-^) ^ ^ S +^ S 2 n + 0(v^). (3.3) 
Inserting such results in (2.30), we get 

e n ~gs(n + ±) + \(3s 2 n 2 + O(y/o) (3.4) 

where 

2~ + \(3 2 + 0(a) s~l + 0(a). (3.5) 
In the limit ct — > 0, we therefore obtain 

e n (/3) = (n + i) + |/3 (n 2 + n + |) (3.6) 

in agreement with [8, 9]. We conclude that in such a limit, the exponential spectrum of the 
general case reduces to a quadratic one. 

It is worth noting that the spectrum corresponding to a = may also be derived directly 
from SUSYQM and shape invariance without resorting to a limiting procedure. Going back 
to the factorization conditions (2.7) - (2.9) and setting a = therein, we are only left with 
a single parameter g = \f3 + \Jl + \f3 2 since s — 1. Among the shape invariance conditions 
(2.15) - (2.17), only the first and third ones survive, namely 

9i+i(9i+i -P) =9i(9i + P) e i+1 = ±(g i + g i+1 ). (3.7) 

The solution of the former is gi + \ = gi + (3, thus giving g% = g + i/3, while the latter directly 
leads to the energy spectrum (3.6) by using the relation e n (/3) = Sr=o e «- 
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This shows that the harmonic oscillator with only a nonzero minimal uncertainty in the 
position is shape invariant under translation of the parameter g. The Hamiltonians hi of 
the SUSYQM hierarchy now reduce to 



h- = 1 



1 + i 2 {3 2 + 2i0yjl + \(3^ P 2 + X 2 + % (i(3 + 2^1 + p 



i = 0,l,2,.... 

(3.8) 

From (2.38) and (2.39), it follows that in this special case rriiU 2 = moo 2 , so that only the 
kinetic energy term in the Hamiltonian Hi is changed. 

3.2 Case a. = (5 / 

There also occurs a simplification in the general formula (2.30) for the energy spectrum 
whenever the dimensionless parameters a and (3 are equal, which means that the original 
parameters a and (3 are related through a = m 2 uj 2 (3. In such a case, we find 

i 1 1 + ct 

7 = k = 1 # = s= 7= 9 = i (3-9) 

Vl — a 1 — a 

and thus 

M = 2# u = t = 0. (3.10) 
After some simple transformation, equation (2.30) can be written as 

e n (q) = \{q + 1){(<7 + l)[n], + 1} = i(g + l)([n], + [n + 1],). (3.11) 

On the other hand, the Hamiltonians hi of the SUSYQM hierarchy now become 

h l = \{q\P 2 + X 2 ) + {q + l)[i\ q }. (3.12) 

From (2.38) and (2.39), it follows that in Hi both the mass and the frequency are scaled 
according to = q~ l m, Ui = q l u. 

The energy spectrum (3.11) is similar to that of the g-deformed harmonic oscillator 

h osc = l(q+l){b,b + } (3.13) 

where b + and b are g-deformed boson creation and annihilation operators satisfying the 
relation [21, 22] 

bb + -qb + b = I. (3.14) 
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Such operators can indeed be used to construct n-g-boson states 

l n >« = ^TrTT l°>« ^ = 0,1,2,... (3.15) 
spanning an orthonormal basis of a deformed Fock space i.e., 

oo 

q (n'\n) q = <J n /, n ^ = i". (3.16) 

n=0 

In (3.15), 1 0) g is the normalized vacuum state, i.e., 6|0) 9 = and ? (0|0) 9 = 1, and the 
g-factorial 

f 1 if n = 

[nU = ^ (3.17) 
1 [nUn-l] s ---[l] s ifn = l,2,... 

is defined in terms of the g-numbers [n] q introduced in (2.31). Since 



b + \n) q = yj[n+l] q \n+l) q b\n) q = J\n\ q \n - l) q (3.18) 

the n-g-boson states \n) q turn out to be the eigenvectors of h osc with eigenvalues e n (q), i.e., 
^osc|n) g = e n (q)\n) q , n = 0, 1, 2, . . . . 

This means that in the a = (3 ^ case, the harmonic oscillator Hamiltonian (2.2) 
with nonzero minimal uncertainties in position and momentum must be reducible to the 
g-deformed oscillator Hamiltonian (3.13). Such an assertion is easily proved by setting 



X = ! v /gTT(6+ + 6) P= | v /gTT(6+ -b) (3.19) 

or, conversely, 

6+ = v!tt (a '- iP) "=vrTT (A ' +ii5) - < 3 - 20 > 

It can indeed be shown from the commutation relation (2.1) that b + and b, as defined in 
(3.20), fulfil the g-commutation relation (3.14) and are such that h = h osc . 

4 Harmonic oscillator eigenvectors in the general case 

To construct the eigenvectors of the harmonic oscillator Hamiltonian (2.2) in the case where 
neither a nor (5 vanishes, one has to resort to a generalized Fock space representation [4, 
5, 7], wherein 

X = y i (q+l)(b + + b) P=^)J^(b + -b) (4.1) 
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are represented in terms of creation and annihilation operators 

b + = - 1 | — X-i^p) b= } (— X + i^P) (4.2) 

satisfying the g-commutation relation (3.14). It is worth noting that in the special case 
where a = (3, considered in section 3.2, the operators (4.2) reduce to the operators (3.20) as 
a consequence of equation (3.9). In the general case where a ^ ft, however, the Hamiltonians 
h and h osc , defined in (2.2) and (3.13), respectively, do not coincide any more since we have 
instead the relation 

h = \(q + 1) (( 7 - [(& + ) 2 + & 2 ] + (7 + {b, & + }) • (4.3) 

Hence the eigenvectors of h are some linear combinations of the n-g-boson states defined in 
(3.15). 

SUSYQM and shape invariance provide us with some prescriptions to construct such 
linear combinations [11, 12, 14]. The ground state is indeed the normalized state annihilated 
by the operator B~(g, s), defined in (2.5), which we shall rewrite here as B~(q,t), i.e., 

B-( q ,t)\MQ,t)) = (Mq,t)\Mq,t)) = 1 (4.4) 

while the normalized excited states can be determined recursively through the equations 

\^n+i(q, t)) = [e n+1 (q, t) - e (q, t)]- 1/2 5 + (g, t) \^ n (q, t,)) n = 0, 1, 2, . . . (4.5) 

where, according to (2.28), t\ = t/q. 

On using (2.21), (2.28), (2.29), and (4.1), the operators B ± (q,t) can be written as linear 
combinations of b + and b, 



B+(q,t) = ^ q -±±{ub + -vb) = ±=K{b+-tb) (4.6) 

B-(q,t) = ] J^(ub-vb+) = -LK(b-tb+) (4.7) 

where K is given in (2.34). 

To work out the explicit form of \tp n (q,t)), n — 0, 1, 2, . . . , it is convenient to use 
the (g-deformed) Bargmann representation of the g-boson operators b + , b, associated with 
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the corresponding g-deformed coherent states [22]. In such a representation the n-g-boson 
states \n) q are represented by the functions 

<Pn(q;0 = -^= (4-8) 

so that any vector \ip) q = Y%?=o c n{q)\n)q £ J~q is realized by the entire function i/)(q;£) = 
J2^=o c n(q)<Pn(q] £)> belonging to a g-deformed Bargmann space B q , whose scalar product 
has been given in [22]. The operators b + and b become the operator of multiplication by a 
complex number £ and the g-differential operator T> q) defined by 

VM«Q (4.9) 

respectively. Hence, B ± (g, t) are represented by 

B + (q,t) = ±=K{S-W q ) B-(q,t) = l=K(V q -tO. (4.10) 

The first equation in (4.4) can therefore be rewritten as 

(v q -tOMq,t;0 = o. (4.ii) 

As detailed in the appendix, the solution of this first-order g-difference equation can be 
easily obtained as 

MQM)=M>(Q,t)Ef (j^ 2 ) (4-12) 
where the g-exponential E q (£) is defined by [23] 

oo tn 
n=0 V L \q- 

and Ao(g, t) is some normalization coefficient, determined by the second condition in equa- 
tion (4.4). 

From (4.8) and (4.13), it follows that the ground state Bargmann wave function can 
also be written as 

00 (Ylv - 11 "\ 1/2 
Mq, t; = A/- (g, t) £ [ [ [2u] J]" ) ^ 2 ^ ; ® (414) 

where 

f 1 if v = 

[2v-lU\={ (4.15) 
\ [21/ - l] q [2u - 3} q ■ ■ ■ [l] q if ^ = 1,2, ... 
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and 

f 1 if v = 

[2uU\ = \ . (4.16) 

\ [2u] q [2u-2] q ...[2] q if 1/ = 1,2,... 

The ground state is therefore a superposition of all even-n-g-boson states. 

The ort honor mality of the functions <fi n {q',£,) in Bq implies that the normalization coef- 
ficient in (4.12) and (4.14) is given by 



oo ro,. ii ii \ 1 / 2 



= [E^gF^j = [i0o(g;-;gV 2 r 1/2 . (4.17) 

Here r s (ai, 0,2, ■ ■ ■ , o-r] &2, • • • , b s ; q, z) denotes a basic hypergeometric function (see [24] 
and the appendix). It can be easily checked that the series in (4.17) is convergent for all 
allowed t 2 values, hence showing that the ground state eigenvector is normalizable as it 
should be. 

Considering next the excited states, we find that in Bargmann representation, the re- 
cursion relation (4.5) becomes the equation 

^n+i(g,t;0 = |[« + i] 9 (i-^r)} ' (£-«Wn(</,*i;0 (^) 

where use has been made of equations (2.33) and (4.10). As shown in the appendix, the 
solution of this equation is given by 

Mq,t;Z)=Mn(q,t)P n (q,t;S)E q 2 L— g 2 ^ . (4.19) 

Here P n (q, t;£) denotes an nth-degree polynomial in £, satisfying the relation 

Pn+i(q, t; = iP n (q, ^ - i^Pn (q, J qtj - tV q P n (q, ^ (4.20) 

with Po(q, t;£) = 1, and Af n (q, t) is a normalization coefficient fulfilling the recursion relation 



K + i(q, t) = I [n + l] q (l - -Lj j I Af n (q, *-j (4.21) 

with J\f (q,t) given in (4.17). 

The solution of equation (4.21) is easily found to be given by 

K(q,t) = {M,!(9- 2n+1 t 2 ;g)ni0o(g;-;g 2 ,g- 2n t 2 )}' 1/2 (4.22) 
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where the symbol (a; q) n is defined in equation (A. 5). 

Solving equation (4.20) for the polynomials P n (q, t; £) looks however more involved. First 
of all, let us remark that this equation may be considered as some g-difference analogue of 
a differential equation satisfied by Hermite polynomials 



H, 



n+1 



2xH n (x) - H' n (x) 



(4.23) 



which can be obtained by combining equations (22.7.13) and (22.8.7) of [25]. As a matter 
of fact, in the formal limit where q — > 1 while t 2 remains different from one^, equation (4.20) 
can be rewritten as 

Qn+i(f; = (1 - t 2 )ZQ n (t; - tQ' n (t; (4.24) 

with Q n (t;£) = P n (l,i;£), because the g-differential operator V q goes to the ordinary 
differential operator d/d£. From (4.23), the solution of (4.24) is found to be 



Q n (t;£) = c n (t)H n [a(m a(t) = c(t) = \i\t{l-V). 



(4.25) 



2t w V 2 

We may therefore regard P n (q,t;^) as some two-parameter deformation of Hermite 
polynomials. As far as we know, such a deformation has not been considered elsewhere. 
In the remainder of this section, we shall therefore devote ourselves to deriving an explicit 
solution to equation (4.20). 

It is straightforward to show that for the first few n values, the polynomials P n (q, t;£) 
are given by 



Pi(q,t;0 = 
P2(q,t;0 = 
Ps(q,t;0 = 



q- 



q 



e-t 



r 



i-^ i-^ 



t 



l--)e 3 --(l + g + rt 

q q 



q J \ q 

For general n values, it is clear from the structure of equation (4.20) that 

Pn(q,t;-£) = (-l) n P n (qM) 



(4.26) 
(4.27) 
(4.28) 



(4.29) 



^By formal limit, we mean a limit that does not correspond to any physical values of the parameters a, 
(3. When q — > 1, we indeed obtain from (2.24) that either a — > or j3 — > 0, which, from (2.11), (2.19), and 
(2.28), implies that either t — > — 1 or t — > 1. 
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which actually agrees with the examples shown in (4.26) - (4.28). We may therefore look 
for a solution of the type 

P n (q,t;0 = jr \[\ + (-l) B - m ]/ B , m (?,f)f n (4.30) 

m=0 

where f n ,m(q,t), m = n, n — 2, . . . , 0(1), are some yet undetermined coefficients. 

Inserting (4.30) in (4.20) and using the property V q £ m = [m] q ^ m ~ 1 , we get the recursion 
relations 

f n+1 , m {q, t) = (l- gra * m+2 ) /n,m-l (q, ^ " *[™ + 1], fn,rn+l (q, ^ 

m = ra-l,n-3,...,0(l) (4.31) 

fn+l,n + l(q,t) = ( X " f ) («' ^) ( 432 ) 

with fofl(q,t) = 1. In (4.31), we have assumed that f n -i(q,t) = for odd values of n. 
The solution of equation (4.32) is given by 

/»,»(?, = ff f 1 - -£r) = f^r; ^ ( 4 - 33 ) 



fe=0 



n 



with ]lfe=o = 1 an d (a; <? 2 )n defined as in (A. 5). As proved in the appendix, the solution of 
equation (4.31) is provided by 

[n] q \ ( t \( n ~ m)/2 ( t 2 2 \ 

fn > m{q > t] = [m] q \[n-m] q \\ {-^m) q ) (n+m)/2 (434) 

where m = n — 2, n — 4, . . . , 0(1). Comparison between equations (4.33) and (4.34) shows 
that the latter can be extended to all allowed m values, i.e., m = n, n — 2, . . . , 0(1). Note 
that the property f n -i(q,t) = 0, for odd n values, can be retrieved by using the usual 
assumption [n] q \ — > oo if n — > — 1. 

Equations (4.30) and (4.34) therefore provide us with the general solution to equation 
(4.20), corresponding to Po(q, t;£) = 1, and they include equations (4.26) - (4.28) as special 
cases. 

To relate the polynomials P n (q, £;£) with some known basic polynomials, let us first 
express them as basic hypergeometric functions [24] . Distinguishing between even and odd 
n values and using equations (A. 4) - (A. 7), we obtain 

ft 2 \ / t V 
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X 201 



t 2 f 2 
-2v 1 2 2v ( i\S 

,^ rT ;g;g ,-q (g-l)y 



(4.35) 



t 



q 4u+1 " y„ +] Vr(s-i). 



X 201 



9' 



-2i/ 



T ;g 3 ;g 2 ,-g 2 ^(g-l)^ 



- - - ^ / (4-36) 

where z/ = 0, 1, 2, . . . . Since in (4.35) and (4.36), the parameter ai of the basic hypergeo- 
metric series is a\ = q~ 2u , where 2v is a nonnegative integer, we deal here with terminating 
scries. 

From equations (4.35), (4.36), and the definition of little g-Jacobi polynomials [24] 

p n (x; a, b; q) = 2 0i(<T n , abq n+1 ; aq; q, qx) (4.37) 

it is now obvious that the polynomials P n (q, t; £) can be re-expressed in terms of the latter 

as 



£2 i f 2 

-q (q-i)—;z>^> ( i 



i ' 1 A 

t q q 4 



(4.38) 



P 2v+1 (q,t;Z) = (q 3 ;q\ 



■q 2 



x Vv 



Vrte-i). 

t2 f 2 
_2i/-l/„ i \ S . „ 1 2 



-^- 1 (?-i)V;^ 



(4.39) 



t ' y ' g 4,+2 

This completes the determination of the eigenfunctions (4.19) of the harmonic oscillator 
Hamiltonian (2.2) in (g-deformed) Bargmann representation. Combining equations (4.8), 
(4.12), (4.19), and (4.30), we can rewrite them as linear combinations of even- or odd-n-g- 
boson wave functions 



(7=0 

X <P2a(q,£) 



V^+i (<?,£;£) = N 2 v+i{q,t) 

(7=0 



mhi(cr,u) 



a— (j, 



f2u,2„(q,t) 



(4.40) 



2^ T77Z o. 1 iihu+l,2fi+l{Q,t) 



a—f-i 



[2(7-2^!!- 



,2v+l 



x ^ 2f 7+i(g,0 



(4.41) 



according to whether n is even or odd. 
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5 Conclusion 



In the present paper, we have determined in a purely algebraic way both the spectrum 
and the eigenvectors of the harmonic oscillator with nonzero minimal uncertainties in both 
position and momentum by availing ourselves of an extension of SUSYQM and shape in- 
variance powerful techniques to the case of the deformed canonical commutation relation 
(1.1). 

In the present context, shape invariance is related to the scaling of some parameter 
t, depending in a complicated way upon the two deforming parameters a, (3 (or the two 
dimensionless ones a, (3) entering the canonical commutation relation. As occurs in other 
examples involving parameter scaling [17], the oscillator spectrum turns out to be expo- 
nential. The supersymmetric partner Hamiltonians correspond to both different masses 
and frequencies. Such an unusual feature is a direct consequence of the deformation of 
the commutation relation and is distinct from the combined potential and effective-mass 
variations that may be effected in the case of a position-dependent effective mass in the 
context of conventional SUSYQM [20]. 

We have proved that whenever one of the deforming parameters vanishes, e.g., a — > 
and P 7^ 0, our exponential spectrum goes to the quadratic one, previously found by solving 
the deformed Schrodinger differential equation [8, 9]. Such a quadratic spectrum may also 
be derived from SUSYQM and shape invariance connected with parameter translation. 

Furthermore, we have shown that when a = j3^0ora = m 2 uj 2 P ^ 0, the harmonic 
oscillator with nonzero minimal uncertainties in both position and momentum reduces to 
the g-deformed harmonic oscillator corresponding to q > 1 [22] and its eigenvectors therefore 
coincide with the n-g-boson states with n — 0, 1, 2, 

Finally, in the general case where O^a^/^^O, we have constructed the oscilla- 
tor eigenvectors as linear combinations of n-g-boson states by resorting to a (g-deformed) 
Bargmann representation of the latter and to g-differential calculus. The ground state can 
be expressed in terms of the g-exponential of an operator proportional to the square of the 
g-boson creation operator, acting on the vacuum, while the excited states contain as extra 
factors nth-degree polynomials in the creation operator. The latter are some two-parameter 
deformations of Hermite polynomials and can also be related to little q- Jacobi polynomials. 
It is worth noting that operators similar, but not identical, to that occurring in the ground 
state eigenvector are familiar in other contexts, such as those of squeezed states [26] and of 
Bose-Einstein condensates [27]. 
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Appendix. Some results used in the determination of 
the harmonic oscillator eigenvectors 

The g-exponential used in the present paper is defined by 



oo tn 
n=0 V L \q- 

where the g-factorial [n] q \ is given in (3.17). For q > 1 (see equation (2.24)), it converges 
for all finite values of £ G C. It is the solution of the g-difference equation 

V q E q (a£) = aE q (a£) a EC (A.2) 

subject to the condition that E q (0) = 1. 

It follows from definition (4.9) and property (A.2) that 

7,p^2x E q ^e)-E q ,{ae) , E#{atfe) - E#{ae) 

= a(g + l)^ 2 « 2 ). (A.3) 

Hence the function ipo(q,t;£), defined in (4.12), satisfies the difference equation (4.11). 

The basic hypergeometric function r s (ai, a 2 , . . . , a r ; bi, b 2 , . . . , b s ; q, z), generalizing the 
conventional one r F s (a\, a 2 , ■ ■ ■ , a r ; b±, b 2 , ■ ■ ■ , b s ; z), is defined by [24] 

r s (ai, a 2 , . . . , a r ; bi,b 2 ,..., b s ; q, z) 

\p ( Q li g)n(Q2; q)n [ [ ' \flr\ QJn r/ \ n ra(n-l)/2l 1+s-r n / * A 

^o(^9)n(&i;g)n(6 2 ;g),---(& s ;g)n K J g J 1 ' J 

where z£C and 

f 1 if n = 

(a; o) n = < . (A.5) 

\ (l-a)(l-aq)---(l-aq n - 1 ) ifn = l,2,... 

In (A. 4), it is assumed that the parameters b±, b 2 , . . . , b s are such that the denominator 
factors in the terms of the series are never zero. Since (q~ m ; q) n — if n — m+1, m+2, . . . , 
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an r (f) s series terminates if one of its numerator parameters is of the form q~ m with m — 0, 
1, 2, .... When dealing with nonterminating series, it is assumed that the parameters and 
the variable are such that the series converges absolutely. 

Some useful relations connecting (a; q) n with the g-double factorials defined in (4.15) 
and (4.16), as well as the latter with g-factorials, are 

\2u] ll= ( g2;g2 )" \2u-l] M = (gA y^M. (A 6) 

L^-- ^ i J«" (g + 1)-- lA - DJ 

We may also note the interesting relations [24] 

(a; q) n+k = (a; g)„(ag ra ; g) fc {a^q 1 ^ q) n = (a; ^(-a" 1 ) V n(n_1)/2 - (A.7) 

Let us now prove that the functions ip n {q,t;£), defined in (4.19), solve equation (4.18) 
provided conditions (4.20) and (4.21) are satisfied. For such a purpose, we shall apply the 
well-known rule for g-derivating a product of functions [23] 

v q f(0g(0 = f(qOv q g(0 + <K0*V(0- ( A - 8 ) 

Substituting equation (4.19) into the right-hand side of (4.18) and using (A. 8) and (A. 3) 
successively, we obtain 



*1 z-2 



(q + l)q n 



^n+i(g,*;0 

= + A4(g,ti)|^n(g,ti;0^ 

-tP n (q,t i; qOV q E q2 (j^^e) -tE? (j^^^ V q P n {q M ^)) 
= |[n + l]^l-^| ' M n {qM)\iPn{q^i)-tP n {qM-,qi) l -^i 

-W q P n {qM^)E q ^[j^^ n e) (A.9) 

which should coincide with equation (4.19) with n + 1 substituted for n. Comparison 
between the right-hand sides of both equations directly leads to (4.20) and (4.21), which 
completes the proof. 

Let us finally consider the solution to the recursion relation (4.31). Substituting equation 
(4.34), where use is made of definition (A. 5), into the right-hand side of (4.31), we get 

fn+i,m{q, t) 
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q n-m+2 J [ m _ [ n _ m _|_ ^ ? (n+m-l)/2 



X 11 1 „2fc+3 

fc=(n-m+l)/2 V y 



nl„! / t \( n_m_1)/2 Vr 1 / t 2 



P^^H + ^n-m + ll.) 



fc=(n-m-l)/2 
(n-m+l)/2 



g2fc + 3 



g(ra+m-l)/2 



n-1 / f 2 



X 



n i 



\ n 2k+3 / 

fc=(ra-m-l)/2 \ y / 

I \ I : <A (A.10) 



H 9 ![n-m+l],!! V q^~^ J \<f n+1 ) (n+m+ i)/2 

where, in the last step, we employ the relation 

[m] q + q m [n - m + l] q = [n + l] q (A.ll) 

and equation (A. 5) again. The final result in (A. 10) coincides with equation (4.34) for 
n replaced by n + 1. We conclude that equation (4.34) provides us with the solution to 
equation (4.31) corresponding to fo,o(q,t) = 1. 
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